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Abstract 

After a brief review of the role three-nucleon forces play in the few-nucleon 
systems, the chiral-perturbation-theory approach to these forces is discussed. 
Construction of the (nominal) leading- and subleading-order Born terms and 
pion-rescattering graphs contributing to two-pion-exchange three-nucleon forces 
is reviewed, and comparisons are made of the types of such forces that are used 
today. It is demonstrated that the short-range c-term of the Tucson-Melbourne 
force is unnatural in terms of power counting and should be dropped. The class 
of two-pion-exchange three-nucleon forces then becomes rather uniform. 
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Introduction 



Three- nucleoli (3N) forces have come under increasing scrutiny recently PJ. Al- 
though these forces are rather weak, they are playing an important role in the theory of 
few-nucleon systems, where computational advances permit calculation of new observ- 
ables that are challenged by experiment [0, [|. The most-recent (second-generation) 
nucleon-nucleon (NN) potentials @, [5| fit the entire NN data base rather well (rival- 
ing phenomenological partial-wave analyses in the best cases) and lead to predictions 
for most 3N observables that are in good agreement with experiment. In a few cases, 
such as the A y puzzle |T|, |2| and the binding energies of few-nucleon ground states 
||, there are inadequacies with this methodology that have focused attention on 
three-nucleon forces. 

All realistic NN forces underbind the triton J7j, and small differences among 
them can be traced to nonlocalities. Three-nucleon forces are incorporated into the 
Hamiltonian and adjusted to achieve the correct triton binding. With this addition 
4 He is properly bound ||, while the two 5 He p-levels have a splitting roughly 30% 
too small 0. Binding of A = 6-8 ground and low-lying excited states is too low ||. 

The best-studied of these problems is the A y puzzle. The calculated asymmetry 
(A y ) in neutron-deuteron and proton-deuteron scattering at low energies is 25-30% 
too small, which looks suspiciously similar to the 5 He problem, since A y is most 
sensitive to spin-orbit forces. A recent analysis of the former problem concludes |I| 
that reasonable changes in the NN force will not resolve the puzzle and that one 
should implement refined 3N force models. Although credible examples of these 



models first began to appear 40 years ago [10], technical problems associated with 
nuclear- force construction hampered the effort, and general acceptance of such forces 
was delayed until it was demonstrated that good NN forces could not reproduce the 
triton binding energy. 

Construction of potentials always involves theoretical choices, since a potential 
is a subamplitude (an off-shell part of an amplitude) that when iterated (in the 
Schrodinger equation, for example) produces observables (on-shell amplitudes or en- 
ergies). The off-shell question has always been a murky one, since it is usually ill de- 
fined. Nevertheless, the same Lagrangian (i.e., the same theory) can lead to different 
potentials, although they should individually produce identical observables. Coupled 
to this is the worse problem of unraveling the underlying strong-interaction physics 
(i.e., deciding on a Lagrangian or equivalent formalism to use). In the early days a 



frequently asked question |l Ij was: how does one account for the off-shell nature of 
(virtual) pions exchanged between nucleons? Faced with such daunting theoretical 
obstacles, all models were simplified. Nonlocality (nucleon-momentum dependence) 
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was typically ignored, for example. The early history of the field is well reviewed in 
Refs. [O, O. 



Since these early beginnings a new formalism |14|, [15], [TJ| [T7[ has been developed 
for implementing strong-interaction physics in low-momentum (for nucleons) regimes: 
chiral perturbation theory (CPT). This technique implements (approximate) chiral 
symmetry (manifested by the quarks in QCD) in constructing the strong-interaction 
building blocks, which are then assembled in all possible ways in the most general 
Lagrangian consistent with the symmetry. At the same time, the entire framework is 
organized with a power-counting scheme. A successful perturbation theory must guar- 
antee that succeeding orders diminish, and chiral symmetry provides the constraints 
mandating that more complex calculations (loops, etc.) should yield progressively 
smaller results, even though strong-interaction coupling constants are not small. This 
scheme also provides a testing mechanism for nuclear interactions: naturalness and 
naive dimensional power counting 



Chiral perturbation theory simplifies the old-fashioned nuclear-physics approach 
of incorporating into a field theory all known meson and baryon resonances with 
energies less than some large (arbitrary) cutoff. All such heavy resonances (with the 
possible exception of the low-lying A isobar, which is ignored here for simplicity) are 
subsumed in short-range (point-like) vertices. In the usual SU(2) approach this means 
that only pion and nucleon fields contribute explicitly, although the entire zoo of heavy 
elementary particles contributes implicitly to the phenomenological constants of the 
theory. 

Two scales that set the strength of the Lagrangian building blocks are f„ ~ 
93 MeV (the pion-decay constant) and A ~ 1 GeV (the large-mass QCD scale). 
Overall powers of A must be negative (i.e., A _A , with A > 0), since they arise from 
the frozen propagation of the heavy states, and interactions in the Lagrangian are 
organized by these powers: C^ A \ Dimensionful coupling constants in this scheme can 
be written as powers of f n and A times dimensionless coupling constants ~ ±1. The 
latter requirement is called naturalness. "Unnatural" implies very small or very large 
(compared to 1) and of either sign. We will use this test later. 



We wish to examine and compare the two-pion-exchange three-nucleon forces 
(3iVFs) that incorporate at least minimal phenomenology from ir-N scattering. There 
are basically four types (plus variants of each that we will not treat): (1) Tucson- 
Melbourne force ]TI[ (the first of this class), based on current-algebra arguments; 
(2) Brazilian force [TJ5[], based on a chiral Lagrangian and a supplemental current- 
algebra constraint; (3) Texas [20 force, based on chiral perturbation theory; (4) 
Ruhr (Pot) force ||21|| , based on non-chiral Lagrangians. Each contains a a-term (or 
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functional equivalent) for s-wave, isospin-symmetric pions, as well as p-wave pions in 
both isospin-symmetric and -antisymmetric configurations, such as might arise from 
virtual A-isobar excitation. We note that the august Fujita-Miyazawa flO| 3iVF con- 
tained equivalents of all these elements (although the s-wave part was dropped) and 
the Urbana-Argonne [22 model contains a conventional Fujita-Miyazawa A- mediated 
force plus an intermediate-range isospin- and spin-independent component. 



Chiral Perturbation Theory 



We will make our comparisons using the framework of CPT, which allows us to 
define the theory in a consistent and transparent way. The relevant parts of the 
leading-order Lagrangian (corresponding to A = 0), are given by [^(], [23| 
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whose three terms correspond to free pions, the free-nucleon energy and Weinberg- 
Tomozawa two-pion interaction, and the usual pion-nucleon interaction. We have 
simplified the nonlinear realizations of the SO (4) symmetry |2U| and dropped terms 
that would have added even numbers of pion fields to all terms with pion fields; we 
do not require such terms in what follows. In addition, the A = 1 Lagrangian, 
is given by p(| 
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where terms with additional pion fields have been dropped, and we have not listed p0| 
three separate spin- and isospin-dependent short-range 3ATF terms (~ (N^N) 3 ) with 
coefficients, e^. We have also ignored isospin violation. Where appropriate we have 
adopted the notation of Ref. |17| and have explicitly incorporated higher-order terms 
resulting from a nonrelativistic reduction of the pseudovector-coupling Born term. 
The phenomenological coefficients q and di must be determined from experiment. 



We have not written down explicit A-isobar contributions above. They are im- 
plicitly included in the phenomenological coefficients. This hides the fact that those 
coefficients that contain tree-level A contributions are expected to be larger than ones 
that do not by a A/(ttt.a — tun) factor. The alternative is to include a A field and 
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count it as a nucleoli field |pD"|. This shifts the nominal order of the isobar effects but 
of course not their numerical value, and it unnecessarily complicates the following 
discussion. 



For later use we also list infinitesimal generators for the (approximate) axial sym- 
metry present in this Lagrangian, where again we ignore terms with more than two 
pion fields: 

7r -> 7T - f n e , (3a) 
N ^N~ i e ' TX7T N, (36) 

J IT 

where e is (a constant) infinitesimal. Under this transformation the three terms 
in Eq. (1) are separately invariant in the limit of vanishing pion mass, as are the 
first bracketed term and each remaining term in £W (j n the same limit). Thus, the 
Lagrangian in Eqs. (1) and (2) is term-by-term (as we have written them) invariant, 
except for the pion- mass and ci-term (also conventionally known as the cr-term): 
—Ami c i — °- 

It is important to note that the Lagrangians are not unique. Redefinition of 
the (unphysical) fields leads to other forms. The form we have chosen satisfies chiral 
constraints in a term-by-term fashion, rather than relying on cancellations between 
sets of terms. It is only important that the chosen form have sufficient generality 
(i.e., enough linearly independent terms). Different forms will then be physically 
equivalent on shell, but will in general be different off shell. Off-shell differences do 



not affect physical processes. Note that the Lagrangian of Ref. ||24|| , which is based 
on a non-relativistic reduction of the relativistic pseudo- vector pion-nucleon coupling, 
used an off-shell extension specified by a continuous parameter, «. Only the choice 
H = 1 corresponds to Eq. (2) and only that choice satisfies term-by-term chiral 
symmetry. Amplitudes calculated using various values of /i correspond to a unitary 
transformation of the Hamiltonian and therefore do not alter physical amplitudes 
(although they are different off shell). We note[F23[] that many of the older papers in 
the field have implicitly adopted different values of \x [viz., -1,0,1]. 

In order to determine the 3ATF to (nominal) subleading order, we need to calculate 
the diagrams of Fig. (1). The two interaction terms in together with the first two 
terms in D- 1 ' are usually called relativistic Born terms, and are separately calculated 
using (the many orderings of) Figs, (la) and (lc), and then subtracting the iteration 
of the one-pion-exchange potential (OPEP) given in (lb). In the static (leading- 
order) limit (mjy — > oo) they have long been known to vanish If one works 
to subleading order one is faced with choices, because different off-shell choices for 
the subtracted OPEP lead to different forms for the 3ATF. Thus, the choice of form 
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(a) 



(b) 



(c) 





(d) 



(e) 



Figure 1: Various three- nucleon- force components that arise in sub leading order in 
chiral perturbation theory, as discussed in the text. 



for OPEP (to order (v/c) ) determines the form of this (Born-term) part of the 3NF. 



The reader is referred to Refs. |2(J and f24fl , where different off-shell choices are 
made. The complete (/i, v) off-shell ambiguity is discussed in the latter reference and 
approximate Lorentz invariance is demonstrated. The former ambiguity arises from 
a nucleon-field transformation (a "chiral rotation") that breaks term- by-term chiral 
invariance, as we discussed below Eq. (2). Different values of [i have been implicitly 
assumed in the past by differing treatments of the Born terms (see Appendix of 
Ref. ||24|| ). The //-dependence arises through differing treatments of the difference 
between (four- vector) q 2 and q 2 (see Eq. (4b) below), and is sometimes called the 
quasipotential parameter. Different quasipotential equations correspond (in part) to 



different values of u, and the values [0,1/2,1] have been commonly used [24]. Different 
values of v correspond to different off-shell amplitudes, but unitarily-equivalent on- 
shell values. Other calculations have ignored part or all of the subleading-order Born- 
term contributions. We will ignore the Born terms in what follows. 

The remaining 9 terms of C^ 1 ' [labeled by Ci,di, e,] generate 3iVFs of the type in 
Fig. (lc) [ci,c 3 ,c 4 ], Fig. (Id) [di,d 2 ], and Fig. (le) [e 1 ,e 2 ,e 3 \. The -k 2 term in £ (1) 
generates contributions of A = 3 size (each time derivative is the same as a nucleon- 
energy difference) and can be neglected. A wide range of physics is subsumed in each 
category. The C3 and C4 terms receive important contributions from A-isobars at the 
blob of Fig. (lc), while a heavy scalar-isoscalar meson would likewise contribute to c\. 
We note that all of the models we will compare contain this important physics, either 
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through phenomenological input or via explicit heavy-particle intermediate states. 



We summarize by noting that the Born term from O- \ the Cj 7r-rescattering 
terms, the di one-pion-exchange terms, and the (purely) short-range e, terms are all 
nominally the same size, although large A-isobar contributions can be expected to 
make some of the terms larger than others. We will not discuss the di and e; terms 
further. This force was first derived in Ref. 120 



Comparisons 



To facilitate comparisons we adopt the familiar framework of the Tucson- Melbourne 
collaboration |lTJ for the Born-subtracted amplitudes 
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a/3 



a + bq-q' + c{q 2 + q' 2 )1 - d (r^ 7 a 3 ■ q x q ') , (4c) 



where 5-functions, phase-space factors, etc., have been ignored, and the invariant 
amplitudes of O, EH| have been expanded in l/m N . 




Figure 2: Contribution to the three-nucleon force arising from a pion emitted by 
nucleon 1 and scattering from nucleon 3 before being absorbed by nucleon 2. 

Equation (4b) is illustrated in Fig. (2), showing nucleon (3) scattering a pion 
emitted by nucleon (1) and absorbed by nucleon (2). The T-matrix for ir-N scattering 
(alone) is denoted and is usually rewritten in terms of F a/S , where a and f3 are 
the isospin labels of the initial and final pions. The final expression in Eq. (4c) 
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holds for pions that have a low momentum (< m n ). Summing over the symmetric 
permutations of (1,2,3) in Fig. (2) leads to the complete three-nucleon potential. 



One easily finds from Eq. (2): 
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J TV 
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where u and u/ are the initial and final pion energies. We have dropped Born term 
contributions to c 2 and C4 in accordance with our earlier discussion. Eq. (5) together 
with the dropped pieces generates the CPT n-N amplitude to 0(Q 2 ). Calculations 
including loops and new parameters at 0(Q 3 ) have also been performed. They have 
been used with different pieces of ir-N scattering data to determine the coefficients 
q. In table (1) we list some of these determinations. Earlier fits |28|, |17], ^9], [3~0| 
were made to different sets of threshold and sub-threshold parameters obtained from 



dispersion analyses of older data. Newer fits |3T[ were made to different phase-shift 
analyses (PSAs), the last two in Table (1) including the more modern meson-factory 
data. The 0(Q 3 ) determinations are consistent with each other when their error 
bars (not shown) are considered, except for ci, which reflects the higher value for the 
cr-term in the newer PSAs. Note that the coefficients c 2 , c 3 , and C4, which receive 
contributions from the A at tree level, are larger than ci, as expected [ j32[| . 



Table 1: Low-energy CPT coefficients in GeV 1 from several recent fits. 



Fit 



0(Q 2 ) 
0(Q 3 ) 
0(Q 3 ) 
0(Q 3 ) 
0(Q 3 ) 
0(Q 3 ) 
0(Q 3 ) 



I 

13 

Tl 



Cl 



c 2 



c 3 



-0.64 
-0.87 
-0.93 
-1.06 
-1.27 
-1.47 
-1.53 



1.78 
3.30 
3.34 
3.40 
3.23 
3.21 
3.22 



-3.90 
-5.25 
-5.29 
-5.54 
-5.93 
-6.00 
-6.19 



2.25 
4.12 
3.63 
3.25 
3.44 
3.52 
3.51 



From the definition (4c) of the (a, b, c, d) coefficients we obtain, 
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Note that there is no c-term, and that the a-term is opposite in sign to the TM result, 
although with C3 < and C4 > 0, b and d are negative and agree with the correspond- 
ing TM signs. A similar result was found in the first of the Brazil-force papers ||19|| , 
where a field-theoretic calculation of isobar contributions was performed. The cr-term 
was not calculated using Feynman rules derived consistently from a Lagrangian, but 
inferred from a tt-N amplitude derived elsewhere. In a later paper, a different off-shell 
amplitude (the one used in the TM calculation) was incorporated. Values of the a — d 
coefficients for popular three- nucleon force models are displayed in Table (2). Note 
that a = a' + 2 m\ c = 1.03/m,r for the TM force. 



Table 2: Low-energy pion-nucleon scattering parameters (with Z-graph [Born] terms 
removed) for a variety of 27r-exchange three-nucleon forces We have also defined a' = 
a — 2 m\ c. The quantities a and a' are in units of m^ 1 , while b, c, and d are in units 
of m~ 3 . 



Three- Nucleon Force 


a' b c d 


Fujita-Miyazawa[|TO|j 
Tucson- Melbourne |TT|, 27 
Brazil JTJ, [Z7| 
Urbana- Argonne 22. 
Texas H|, H 
Ruhr(Pot)|2| 


0.0 -1.15 0.0 -0.29 
-1.03 -2.62 1.03 -0.60 
-1.05 -2.29 1.05 -0.77 

0.0 -1.20 0.0 -0.30 
-1.87 -3.82 0.0 -1.12 
-0.51 -1.82 0.0 -0.48 



Given that CPT is a comprehensive approach to calculating strong-interaction 
physics based on chiral symmetry and subsumes current algebra |17|, j33H , how can 
the CPT (corresponding to derivatively-coupled pions) amplitude [Eq. (5)] and TM 
amplitude [Eq. (4c)] differ? 



We answer that question by noticing that the difference resides only in terms that 
vanish when the pions are on-shell (as we shall see in Eq. (11)). We return to the 
earlier off-shell discussion, and follow closely the approach of Ref. p|. Off-shell 
amplitudes are not unique, and in a field-theoretic calculation they depend on the 
fields chosen to represent pions and nucleons. Our form was chosen to satisfy chiral 
symmetry term-by-term, thereby attaining manifest power counting. Current-algebra 
constraints at certain off-mass-shell points [26| are not satisfied by our isospin-even 
Ti-N amplitude, [F a/3 = 5 a ^F^ + ' + •••]. These points all correspond to vanishing 
(four- vector) q ■ q', as well as u and uj' (to the order we work). Consequently, we can 
ignore the C2- and C3-terms in Eq. (5) and concentrate on the remaining term, which 
can be written in the form 



-^cpt — "H > ( 7 ) 

J TT 
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which holds everywhere. 



Again following Ref. | 33| , we redefine the pion field as 

tt' = tt (1 - J- N^N) , (8) 

and work only to order A = 1 (since o ~ 1/A). Substituting Eq. (8) into £(°\ we 
generate the extra terms 



A£« = 



A f N(n' ■ Dtt' + mjy 2 ) - -^-A^ a • V[r • k'N^N] N 

2 



+ 



(9) 

The last term involves four nucleon fields and is not immediately required. The three 
terms in and involving a and two nucleon fields lead to 



F^^^^ + q'i-ml), (10) 

TV<I 7T 



which agrees with Eq. (7) at any on-shell (e.g., Cheng-Dashen) point (q = q — m 2 ), 
but vanishes at the Adler points (q 2 = m 2 ,g' 2 = 0, and q' 2 = m 2 ,q 2 = 0), and at 
the Weinberg point (q 2 = q' 2 = 0) has the value: F$ = —a/f 2 . Equation (10) 
therefore agrees with the usual current- algebra constraints [JIT], |2(|, as does our entire 
amplitude, F^ + \ in the new pion-field basis. Thus, there is no conflict here between 
CPT (with derivatively-coupled pions) and an approach based on current algebra. 
The only difference is in the choice of fields used to specify the chiral Lagrangian, 
and observables calculated for physical processes must be identical. 

In the TM approach it was noted that rewriting Eq. (10) in terms of inverse pion 
propagators, 

F^ = Y2 +dr7^ 2 -™'+?' 2 -^)> ( n ) 

Jn ' n nJiT 

allows cancellation of the pion propagators in Fig. (2). The first (constant) term 
reproduces Eq. (7) (-^cpt)- This rearrangement amounts to undoing the field trans- 
formation in Eq. (8) that led to Eq. (9), and leads to an effective a-term (a') that 
has a common sign for all models: a' = a — 2 m 2 c. Cancelling the inverse propagators 
in the second term in Eq. (11) leads to a new short-range-plus-pion-range 3AF: 

However, a three- nucleon force of the same type is generated by the last term in Eq. 
(9), comprised of four nucleon fields and one pion field, together with the last term 
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in Eq. (1): two graphs as in Fig. (Id) give 




q 2 + rri: 



1 




q ' 2 + rri. 
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7T 



T )T! -T 2 . 



(13) 



This is exactly equal in size and opposite in sign to the new short-range contribution 
from the off-shell extrapolation of the tt-N amplitude, Eq. (12). This cancellation 
is to be expected, since our original (chiral) Lagrangian produced no such terms to 
start with, and we have just been rearranging terms since then. In summary, the 
TM approach used a current-algebra representation of the amplitude, performed an 
implicit field redefinition to our (CPT) choice of fields, which resulted in an extra 
short-range term in their result. Why did they have an extra term and we do not? 

The TM calculation was predicated upon current- algebra constraints on the off- 
shell tc-N scattering amplitude, which can be reproduced in the CPT approach, as 
well, as we have demonstrated. It is not enough, however, to worry about just 
that scattering amplitude, if one constructs a 3iVF. To incorporate all of the chi- 
ral constraints into the three-nucleon force, current-algebra constraints on the pion- 
production amplitude from two nucleons would also be necessary (leading to the last 
term in Eq. (9)): a daunting task in the current- algebra approach of TM, but one 
that is unnecessary in our approach. We emphasize that a detailed analysis of the 
off-shell region of n-N scattering (for example) is equivalent to a particular choice of 
fields, and (while interesting) is not necessary for constructing a 3NF. 

If one uses the pion- field redefinition in the symmetry generators, Eq. (3), one 
finds that the entire Lagrangian maintains its original symmetry, but that AC^ 
generates new non-invariant terms that cancel against additional contributions from 

£(0) 

(via the new term eoN^ N /m^f^ in the pion generator). One might presume that 
since all of these terms violate chiral symmetry this poses no problem. Unfortunately, 
chiral-symmetry-breaking terms must vanish in the chiral limit. The additional terms 
in Eq. (9) (being just a redefinition of fields) exactly cancel each other in any on- 
shell amplitude. Individually, the two terms do not vanish in this limit because the 
presence of the l/m^ in Eq. (9) removes the implicit m 2 in a, and cr/m 2 does not 
vanish in the chiral limit (m n — > 0). Reiterating, the structure of the additional terms 
in Eq. (9) means that they must individually vanish in that limit, or the entire set 
of terms must be kept to allow for exact cancellations between them to restore the 
proper limit. Because the TM approach (implicitly) kept only the first term in Eq. 
(9), that limit could not be guaranteed for the three-nucleon force. Another way of 
saying the same thing is that dimensional power counting (naturalness) is not satisfied 
for the individual terms in Eq. (9). 

One can check this conclusion by dimensional power counting. An interaction of 
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the form of the last term in Eq. (9) is chiral-symmetry breaking; if it alone is to 
be kept, it has to be implicitly proportional to m%, and hence is nominally an 
term. Such £® coefficients have a generic size im^/^A 3 , where the dimensionless 
coefficient x should be of order 1. If we equate this to gA& /2m^f 3 (the coefficient 
of the last term in Eq. (9)), we obtain x ~ gA&A 3 /(2 m\) ~ 100, which is vastly 
unnatural. The unnatural coefficient [g A cr /2rri^ fl\ is entirely the result of ojm^ 
having a finite symmetry limit. 

We recommend that the short-range c-term in the TM force be dropped (but 
the full value of a! in Table (2) retained; note that the proper power counting has 
been maintained in a! by the factor of preceding c in the definition of a' = 
a — 2m^c, where now each term in this definition vanishes in the chiral limit). This 
had been previously advocated by the Brazil group for reasons having nothing to 
do with symmetry. We note that the d\ and d2 terms in the Texas force are also 
short-range in one pair of nucleons and of pion range in the other. These parts of 
that force (and the corresponding terms in the Lagrangian) satisfy chiral constraints, 
as does a fully short-range force of the generic type contained in the UA 3iVF, and 
shown in Fig. (le). 

In summary, we have briefly reviewed the class of "realistic" three-nucleon forces. 
We have demonstrated that the short-range c-term of the TM approach is unnatural 
and should not be kept. 
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